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INTRODUCTION 


Me la l matrix composites have relatively low densities and yet can be very 
strong and stiff. They can be used in more aggressive environments than resin 
matrix composites, and their weight saving potential is well documented. 

However , most fibers are elastic and their composites tend to be very notch 
sensitive. Ballistic impacts of unidirectional boron/aluminum caused damage 
that acted like a sharp notch, causing a large loss of strength [1]. Residual 
strengths were predicted from the fracture toughness (stress intensity factor at 
failure) by treating the impact damage as an equivalent crack. Thus, fracture 
toughness is a very important property for metal matrix composites. 

The fracture toughness of composites depends on fiber and matrix 
properties, fiber orientations, and stacking sequence. There are far too many 
combinations of fiber, matrix, and layup to evaluate experimentally. Thus, some 
analytical guidance is needed to select fiber, matrix, and layup to give maximum 
fracture toughness for a given strength and stiffness. Accordingly, a maximum 
strain criterion was used to derive a general fracture toughness parameter that 
is independent of laminate orientation [2]. The fracture toughness can be 
predicted from the parameter using the elastic constants of the laminate and the 
fiber failing strain. These properties are readily obtainable. This method 
gave good results for resin matrix composites [3-5] but not for boron/aluminum 


The boron/aluminum specimens in reference [2] contained central crack- like 
slits and were made with various proportions of 0° and ±45° plies, including 
unidirectional and [±45]2 S laminates. Widespread yielding of the aluminum 

matrix caused the compliance to be very nonlinear, making use of the elastic 
constants to predict fracture toughness and strength erroneous. The degree of 
nonlinearity, and hence the error, increased with the proportion ±45° plies. 
Hence , in reference [2], the stress intensity factor was replaced by a strain 
intensity factor and the singular strain field given by the Theory of Elasticity 
was assumed to be valid, eliminating the elastic compliance from the equation 
for the general fracture toughness parameter. The strain intensity factor was 
derived for a uniaxially loaded specimen with a central crack. Thus, the strain 
intensity factor at failure, rather than the fracture toughness, was predicted. 
Then failing strains were predicted for the various crack lengths and strengths 
were calculated using stress-strain curves. The predicted strain intensity 
factors at failure and strengths were in good agreement with those from the 
experiments, even for the laminate with only ±45° plies. This approach should 
be valid for other metal matrix composites that have continuous fibers. 

More recently, numerous people have calculated fiber stress concentration 
factors at a crack tip in metal matrix composites and accounted for yielding of 
the matrix. The fiber stress concentration factor is the ratio of fiber stress 
at the crack tip to the fiber stress away from the crack. Reedy [6,7] and Goree 
et al [8,9] used a discrete model for unidirectional boron/aluminum composites, 
and Johnson and Bigelow used an elastic-plastic finite element model [10,11] for 
boron/aluminum and silicon-carbide/aluminum laminates with various proportions 
of 0°, ±45°, and 90° plies. Reedy [7] and Goree and Jones [8] also determined 
how changing the aluminum yield strength affected the fiber stress concentration 
factor for unidirectional boron/aluminum. Post et al [12] measured the strains 
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flt the crack tip of a [0/±45]<, boron/aluminum specimen using moire 

interferometry. Goree et al [13] also calculated the fiber stress concentration 
factor and accounted for the effect of off-axis plies (other than 0° plies). 

For both the discrete and continuum modeling [6-11,13], the specimen was 
assumed to fail when the stress in the fiber at the crack tip exceeded its 
tensile strength. It is important to note that the fiber stress failure 
criterion and the maximum strain criterion in reference [2] are equivalent since 
the stress -strain behavior of the fibers is linear. Strain equations were used 
in reference [2] because they are simpler than lamina (fiber) stress equations. 

In this paper, the work in reference [2] is reviewed and related to other 
relevant and more recent work. In the first section of this paper, the material 
and experimental procedure of references [2,14] are described briefly for the 
convenience of the reader. The emphasis in this paper is on the analysis method 
and not the experimental data. The reader is referred to references [2,14] for 
more details on the material, fracture test results, and tensile stress-strain 
behavior. In the second section, the nonlinear stress -strain behavior of 
boron/aluminum is reviewed to illustrate how predictions using elastic constants 
result in overestimations of unnotched strength. More details are given in 
reference [14] for the unnotched stress -strain behavior and mechanical 
properties. In the third and fourth sections, the general fracture toughness 
parameter and the strain intensity factor for a uniaxially loaded specimen 
containing a central crack are reviewed. In the last section, the experimental 
values of the general fracture toughness parameter are presented and the 
predicted values of fracture toughness, strain intensity factor at failure, and 
strength are compared with experimental values. 
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half-length of crack or crack-like slit, m 
factor in Ramberg-Osgood equation, Pa '*' 

Young's modulus, Pa 
Young's modulus of fibers, Pa 

ultimate tensile strength of laminate (unnotched) , Pa 
ultimate tensile strength of fibers, Pa 

stress intenstiy factor at failure or fracture toughness, ?ajm 
elastic stress intensity factor at failure, Pa7m 
strain intensity factor at failure, Jm 
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elastic strain intensity factor at failure, Jm 
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length of specimen, in 

exponent in Ramberg- Osgood equation 

general fracture toughness parameter, Jm 

polar coordinates 

gross laminate (applied) stress in y-direction at failure (strength), 
Pa 

width of specimen, m 

fiber or ply orientation angle (relative to loading axis) 

angle of principal load-carrying plies 

ratio of number of 0° plies to total number of plies 

axial strain 

far- field (remote or applied) strain at failure 

ultimate tensile strain of laminate (specimens without crack-like 
slits) 

ultimate tensile strain of fibers 
Poisson's ratio 

dimensionless material constant 

inherent crack lengths calculated from stress and strain, 
respectively, in 

axial stress, Pa 


Suhsc r i pt s : 

x,y Cartesian coordinates (The y-direction corresponds to the axial 

loading direction of the specimen or laminate.) 

1,2 principal ply coordinates (1 refers to fiber direction) 
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MATERIAL AND EXPERIMENTAL PROCEDURES 


The boron/ aluminum composites in references [2,14] consisted of 0.142-mm- 
diameter (0.0056-in.) boron fibers in 6061 aluminum. The composites were tested 
in the as -fabricated condition. The yield strain of the aluminum was only about 

0. 0005, which corresponds to a strength of 34 MPa (5 ksi) . The laminate 
orientations were [0 ] 6 , [0 2 /±45 ] g , [±45/0 2 ] s> [0/±45] s , and [±45]^. The fiber 

volume fractions were 0.50 for the [0] 6 laminates and 0.45 for the others. The 

mechanical properties and stress-strain behavior are described in reference 
[14]. For the convenience of the reader, pertinent mechanical properties are 
given in Table I. 

The middle cracked fracture specimens of reference [2] are shown in figure 

1. The specimens were 19.1, 50.8, and 101.6 mm (0.75, 2.00, and 4.00 in.) wide 
and at least twice as long as wide. The crack- like slits were cut with an 
electrical-discharge process. The specimens were tested in a hydraulic, servo- 
controlled testing machine. The load was programmed to vary linearly with time 
at a slow rate, typically about 2 minutes to failure. Far-field strains were 
measured with strain gages that were located a distance from the slit of at 
Least three times the slit length. 

The unnotched tensile specimens of reference [14] were rectangular in shape 
and were 19.1 mm (0.75 in.) wide and 254 mm (10 in.) long. The specimens were 
tested in the same testing machine as the fracture specimens and with a similar 
load rate. 


NONLINEAR STRESS -STRAIN BEHAVIOR OF BORON/ALUMINUM 


Composites with fibers all in one direction do not have adequate strength 
and stiffness in the transverse direction for most applications. Thus, most 
composite laminates also contain off-axis plies, plies with fibers at some angle 
« to the loading direction, such as ±45° and 90°. When the yield strength of 
the matrix is exceeded, the matrix carries a lower proportion of the load than 
linear-elastic lamination theory predicts, particularly in the off-axis plies. 
The effect of yielding on the stress-strain behavior of the five boron/aluminum 
laminates is shown in figure 2. The specimens were loaded in the direction of 
the 0° fibers. Yielding of the aluminum causes the curves to be nonlinear. The 
degree of nonlinearity increases with the proportion of ±45° plies. 

For convenience, the stress-strain curves in figure 2 were calculated with 
the following Ramberg-Osgood equation. 
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The values of b^ and n , which are given in Table I, were determined in 
reference [14] by regression analysis of experimental data for 4 to 11 specimens 
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of each laminate orientation. The double subscript notation in equation (1) is 
taken from reference [14] to be consistent. 


do 

The effect of yielding on the tangent modulus is shown in figure 3. 

y 

The tangent modulus , which is divided by the elastic modulae in Table I , 

was calculated from the derivative of equation (1). For small applied stresses, 
the Ramberg-Osgood equation does not model the tangent modulus well for 
laminates that contain 0° plies. The actual yield strengths are larger than 
indicated in figure 3. For example, the yield strength of the [0]^ laminate was 

about 117 MPa (17 ksi) , which corresponds to the aluminum yield strength of 34 
MPa (5 ksi). (The laminate stress is 3.4 times the matrix stress.) The yield 
strength of the [ ±4 5)25 laminates was lower, about equal to that of the 

aluminum. The yield strengths of laminates with both 0° and ±45° plies varied 
between those for [0) 6 and [±45 ] according to the proportion of 0° plies. 

With increasing stress, the tangent modulus in figure 3 decreases and 
asymptotically approaches a value that corresponds to the situation where 0° 
fibers carry all the load. One can infer from the curves that, for small 
stresses in the elastic range, the aluminum matrix of the [0]g laminates carries 

about 20 percent of the load, and the ±45° plies of the other laminates carry 
from 50 to 100 percent of the load, depending on the proportion of ±45° plies. 
Therefore, yielding reduced the load-carrying potential of the laminates 20 to 
nearly 100 percent, depending on the proportion of ±45° plies. 

The average failing strains for the laminates that contain 0° plies are 
plotted in figure 4. The failing strains were essentially equal for the 
different laminate orientations, 0.0076 on the average. For a fiber modulus 

- 400 GPa (58 Msi) , a strain of 0.0076 corresponds to a stress in the 0° fibers 
of 3.04 GPa (441 ksi), which is typical of the strength of these boron fibers. 
Thus, failure of these laminates coincided with failure of the 0° fibers, which 
carried virtually all of the load because of yielding. 

For linear elastic stress-strain behavior, strengths of the laminates with 
0° plies are given by 


e P E 
tuf y 


( 2 ) 


where « ^ " F tuf^f t *' e tens il e failing strain of the fibers and F tu f is 

the tensile strength of the fibers. However, as shown in figure 5, strengths 
calculated with equation (2) and e ^ - 0.0076 are much too large because of 

yieLding, especially for laminates with ±45° plies. 


On the other hand, the unnotched strengths can be predicted using the 
stress-strain curves rather than elastic constants. Assuming that the strains 
in the 0° and ±45° plies are equal at failure, the strength can be written as 


F - 1541A + 143.8(1 - A) 


(3) 
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where A is the proportion of 0° plies and the first and second terms give the 
portion of load carried by the 0° and ±45° plies, respectively. The factors 
] Vi ! and 141.8 MPa were calculated with equation (!) for (lie |()| f and [ | /, * j j ^ ^ 

laminates, respectively, assuming = 0.0076. Equation (3), which is plotted 

in figure 5 as the dashed line, is in good agreement with the data. 

Strength for the [+45 unnotched tensile specimens is not shown in 

figures 4 and 5 because the unnotched tensile specimens failed along a line that 
was 45° to the loading axis, indicating that failure was related to the maximum 
shear stress rather than the maximum tensile strength of the fibers. However, 
the path of the failure in specimens with cracks was mostly through the net 
section, indicating that failure from the crack-like slits was related to the 
maximum tensile strength of the fibers. 


REVIEW OF GENERAL FRACTURE TOUGHNESS PARAMETER * 


Failure of a composite containing a crack is precipitated by failure of the 
principal load-carrying plies, much as in the case of unnotched laminates. The 
principal load-carrying plies are generally the ones with fibers most oriented 
with the applied load (smallest angle a in figure 6), typically 0° plies. At 
failure, the fiber strains ahead of a crack tip (0 = 0 in figure 6) in a 
specially orthotropic laminate under plane stress and mode I conditions [2] can 
be written 


‘lc - Q c (2 " r)1/2 ♦ B 0 + V 1/2 + V 3/2 + ••• (4) 

where r is the distance from a crack tip. The coefficient Q c is given by 

n _5tt£ 

% E y (5) 
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and Kq is the value of the stress intensity factor at failure or the fracture 

toughness. The angle a is the angle that the principal load-carrying fibers 
make with the y-axis in figure 6. For all the laminates containing 0° plies, 

« “ () ; and, for the [±45) 2s laminates, a* - 45°. 


It was assumed that the principal load-carrying plies fail when the fiber 
strains given by equation (4) exceed a critical level. Thus, Q c is a constant 
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at failure. The coefficient was referred to as a general fracture 

toughness parameter [2] because it is independent of laminate orientation. As 
noted previously, this failure criterion is equivalent to a fiber stress 
criterion. It follows from equation (5) that the fracture toughness is 

proportional to the Young's modulus and the nondimens ional factor f. The 

factor £, which is given in Table II and plotted in figure 7, is not strongly 
affected by the proportion of 0° plies for the boron/aluminum laminates. Both 
E and £ can be calculated using lamination theory. It was shown in 


references [3,4] for many resin matrix composites with different laminate 
orientations that the value of Q increased in proportion to the failing 

strain of the fibers e The constant of proportionality was 1.5 7mm on the 


average. Thus, fracture toughness can be predicted with equation (5) from the 
c 1 a s t i c constants and the failing strain of the principal load-carrying fibers. 


A STRAIN INTENSITY FACTOR 


Predicting fracture toughness with equation (5) is analogous to predicting 
unnotched strength with equation (2). In either case, the prediction will be in 
considerable error if the compliance is very nonlinear. On the other hand, a 
strain- intensity- factor approach can be used to predict strengths of cracked 
specimens similar to the strain approach used to derive equation (3) for 
strengths of uncracked specimens. A review of the derivation in reference [2] 
of a strain intensity factor for uniaxially loaded specimens containing a 
central crack follows. 

The fracture toughness or stress intensity factor at failure for a 
homogeneous orthotropic sheet of infinite extent, containing a central crack of 
Length ?a, can be written 


Kq - S^Jn(a 4 p o ) sec (jra/W) (6) 

where S c is the applied stress at failure perpendicular to the crack and 




(7) 


The constant p Q was added to crack length so that S c - F , the unnotched 

strength, when a - 0. The secant function corrects for finite width. It has 
been verified for linear-elastic, orthotropic laminates with 2a/W <0.5 using 
finite element analyses. For example, see reference [15]. 

A more convenient form of equation (6) is 
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( 8 ) 


k q ” V (1 


■--.Qe.. r l/2 
nsFl } 

tu 


who re 


Kq^ - S^Jna sec (jra/W) 

For a uniaxial applied stress, the strength is 

S - e E 
c c y 


(9) 


( 10 ) 


where £ c is the far-field strain at failure. Substituting equation (10) into 

( 6 ), 


K_ ** E K _ 

Q y «Q 

where is the strain intensity factor at failure given by 


(ID 


K €Q - tjn (a + P e ) sec (jra/W) 


- ^> 2 
tu 


(12) 

(13) 


;md € tu is the fai l in 6 strain of the uncracked laminate. For laminates with 

0° plies, f tu ” c tuf’ However . in this paper, all calculations of were 

made using the values of e in Table I. 

tu 

Similar to equation (6), equation (12) can be written 




K _ (1 
eQe 


AqS v-1/2 
2 ' 

7ra e 

tu 


(1A) 


where 


K £Qe ” f c Jna sec (,ra / W) 
Substituting equation (11) into (5) , 


W 


(15) 


(16) 
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Equation (16) relates the level of strain at failure in the principal load- 
carrying plies ahead of a crack tip to the far- field strain. It was assumed 
that the level of strain at a crack tip, when normalized by the far-field 
strain, is not affected by widespread yielding or a nonlinear compliance, which 
is equivalent to assuming that the factor £ is not affected by yielding and 
can be calculated with the elastic constants. 

The strains ahead of a crack-like slit in a [0/+45],, laminate were measured 

by Post et al [14] using moire interferometry. Compared to the change in 
compliance, the level of strain in the 0° fiber direction normalized by the far- 
field strain was relatively constant with increasing applied load for points at 
Least 0.19 mm (0.0076 in.) beyond the crack tip. This distance corresponds to 
the typical spacing of boron fibers, which have a diameter of 0.14 mm (0.0056 
in.). However, for a point at the crack tip, the fiber strain concentration 
factor decreased dramatically with increasing far-field strain, from 8 to 3 . 

In unidirectional laminates, the discrete models [6-9,13] predict a similar 
reduction in fiber stress or strain concentration factor due to yielding of the 
aluminum in shear at the crack tips. Since the fibers have linear-elastic 
behavior, the fiber stress concentration factor is equal to the fiber strain 
concentration factor. The yield zone is less than a fiber diameter in width and 
extends away from a crack tip parallel to the fibers. The length of the yield 
zone at failure can be several times the length of the crack. However, in 
laminates with both 0° and ±45° plies, the ±45° plies bridge the yield’ zone, 
greatly reducing the shear stresses and the length of the yield zone [13]. 
Consequently, yielding in shear at the crack tip reduces the fiber strain 
concentration factor much more for a unidirectional laminate than for a laminate 
with both 0° and ±45° plies, as will be evident in the experimental data. 

It is important to note that the effect of yielding on fracture toughness 
can depend on specimen type. In reference [16], Reedy attributed differences 
between fracture toughness values for center-cracked, three-point-bend, and 
compact specimens made from unidirectional boron/aluminum to widespread 
yielding. Also, the fracture toughness values for three-point-bend specimens 
were different for different sizes, and the compact specimens split at the ends 
of. the crack-like slit and failed without breaking fibers. In contrast, the 
center-cracked and three-point-bend specimens failed along a line coincident 
with the slit. The behavior of these specimens would probably have been less 
deviant had the laminates contained off-axis plies in addition to the 0° plies. 

Also, it is interesting to note that the fiber stress concentration factor 
at the crack tips is reduced far more by splitting than by yielding [9,13,17], 
Thus, the singular strain field given by equation (4) is not valid at all’ when 
long splits develop at the crack tips. Resin matrix composites tend to crack or 
split at the crack tips more so than metal matrix composites. 
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RESULTS AND DISCUSSION 


General Fracture Toughness Parameter 

Values of Q c were calculated for each specimen using equations (14)-(16) 

and plotted against slit length in figures 8(a)-8(e). Average values are 
plotted in figure 9 for each laminate orientation and specimen width. For the 
[0) 6 specimens in figure 8(a), the values of Q c increase with slit length. As 

noted previously, shear yielding at a crack tip reduces the fiber stress 
concentration factor and thus elevates strength and far-field failing strain. 

For a given applied load, the length of the yield zone and hence the reduction 
in the fiber stress concentration factor was shown [8,9] to increase with 
increasing slit length. In fact, for the long slits, the yield zones were long 
enough to alter measurements of remote strain [2], The average Q c values for 

[ 0 ] 6 in figure 9 increase with increasing specimen width because slit length 

increases with specimen width. For a given slit length, the values of Q c in 

figure 8(a) do not appear to increase with increasing specimen width. 

On the other hand, the values of Q for the [0„/±45] [±45/0 o ] , and 

c 2 S 2 S' 

|0 ? /±45]g laminates in figures 8(b) -8(c) do not increase noticeably with slit 

length. As noted previously, the fiber stress concentration factor is reduced 
less by shear yielding at the crack tips when laminates contain both 0° and ±45° 
plies than when laminates contain only 0° plies. However, the average values of 
Q c in figure 9 for these laminates do increase with increasing specimen width, 

and the increase is greater with increasing proportion of ±45° plies. 

For [±45 ] 2s specimens with short slits, the far-field strains at failure 

reported in reference [2] were very large, much as those for the unnotched 
specimens. For wide specimens, the resulting values of Q c were so large that 

a larger scale had to be used for figure 8(e) than for figures 8(a) -8(d). On 
the other hand, far-field strains at failure for the wide [±45 ] 2g specimens with 

long slits were less than one-tenth those with short slits, causing the Q 

values in figure 8(e) to decrease with increasing slit length. In fact, for 
wide specimens with the longest slits, far-field strains at failure for [±45]„ 

specimens were as small as those for specimens containing 0° plies. For this 
reason, the average values of Q £ in figure 9 vary greatly with specimen width. 

For [ ±4 5 ] 2 g specimens with long slits, the values of Q c in figure 8(e) 

are considerably larger for the 101 . 6-mm-wide (4. 00- in.) specimens than for the 
80. 8 -mm- wide (2.00-in.) specimens. Notice how the unnotched strengths in Table 
I increase with specimen width. Therefore, the increase of Q c values with 

specimen width for [0 2 /±45] s> [±45/0 2 ] g , [0 2 /±45] s , and [±45 ] 2g laminates is 

probably related to the ±45° plies. 
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It was previously noted that the unnotched [ ^ ] 2 s s P ec ^- mens failed on a 

plane oriented at 45° to the applied load, indicating that the failure was 
associated with the maximum shear stress. On the other hand, the [±45] 2S 

spec linens with slits failed along a line that was more or less coincident with 
the slit, which is consistent with tensile failure of the fibers. 

As shown in figure 9, an average value of Q c for laminates containing 0° 

plies and for specimens wider than 19.1 mm (0.75 in.) is 0.0121 7mm. It was 
shown [3,4] for many composite materials and laminate orientations that the 
value of increased in proportion to the failing strain of the fibers € tu f* 

The constant of proportionality was 1.5 7mm on the average. Taking € tu f * 

0.00760, and - 0.0121 7mm f° r the boron/aluminum, Q c / € tu f “ 1-59 7mm, which 

is reasonably close to 1.5 7mm. 


Strain and Stress Intensity Factors at Failure 


The values of calculated with equations (14) and (15) were averaged 

for each specimen width and plotted against the percent of 0° plies in figure 
10, as well as recorded in Table II. Values predicted using equation (16) with 
the £ values in Table II and * 0.0121 7mm are plotted as a line for 

comparison. Equation (16) is in good agreement with the test data except for 
the [±45 ] 2s laminates. Again, for the [±45]^^ laminates, the very large far- 


field strains at failure for wide specimens with short slits elevated the 

values. Notice that the predicted value of is relatively constant with 

increasing proportion of 0° plies, indicating that the far- field strains at 
failure for the various laminate orientations are also nearly equal for a given 
slit 1 ength [ 2 ] . 


Average values of fracture toughness or stress intensity factor at failure 
calculated with equations (8) and (9) are plotted similarly in figure 11, as 
well as recorded in Table II. The solid line was predicted using equation (5), 
the elastic constants and £ values in Table II, and — 0.0121 7mm. Except 

for the [0]^ laminates, the solid line is far above the test data, much as in 

figure 5 for the unnotched tensile strength. In effect, widespread yielding 

caused the laminates with and without slits to be weaker than linear- elastic 

theory predicted. The dashed line was also calculated with equation (5) except 

that E was replaced by the secant modulus F^ /e^ . The test data and the 
y tu' tu 

dashed curve are in fairly good agreement. 

For silicon-carbide/aluminum laminates, Johnson and Bigelow [11] predicted 
trends similar to those in figure 11 using an elastic-plastic finite element 
model. For laminates with both 0° and off-axis plies, they found that the fiber 
stress at the c rack tip for a given applied laminate stress or load was more 
when yielding occured than when the material was elastic, indicating that 
widespread yielding causes a cracked laminate to be weaker than predicted by 
linear-elastic theory. However, for a unidirectional laminate, they found that 


11 



the fiber stress at the crack tip was less when yielding occured than when the 
material was elastic, indicating that shear yielding at the crack tips causes a 
cracked laminate to be stronger than predicted by linear-elastic theory. 
Therefore, the analysis in reference [11] and the results of this paper show 
that shear yielding at the crack tip increases the strength of a unidirectional 
laminate, whereas widespread yielding decreases the strength of a laminate with 
0° and off-axis plies. 


Strength Predictions 


Values of were predicted for each laminate orientation using equation 

(16) with Q c - 0.0121 7mm and the £ values in Table II. Then, far-field 

strains at failure were calculated using equations (14) and (15) with Q c - 

0.0121 7mm and the e values in Table I, and strengths were calculated using 

equation (1), the Ramberg-Osgood equation. This approach is analogous to using 
equation (3) to predict unnotched strength. Strengths were also predicted 
directly with the stress intensity factor using equations (5), (8), and (9) with 
Q c ” 0.0121 7 mm , - F tu /e tu > and the values of F fcu and e tu in Table I. 

Both predictions of strength are plotted in figures 12(a) -12(e) along with the 
experimental strengths for all five laminate orientations. The strengths were 
multiplied by the secant correction factor to make the strengths for the 
different specimen widths coalesce. For this reason, only one curve is shown 
for the direct predictions of strength using the stress intensity factor except 
for the [±45]2 <j specimens in figure 12(e). Here, a different unnotched strength 

was used for each specimen width. On the other hand, the strength curves 
predicted from the strains do not exactly coalesce. However, except for the 
[ J.45 ] 2 s specimens, the curves do not differ significantly for 2a < W/2. 


Strengths were not predicted for 2a > W/2 since the accuracy of the secant 
correction factor is questionable. The strengths predicted from the strains and 

the experimental strengths are in fairly good agreement, even for the [±45] 0 

Z j 


specimens. The strengths predicted directly with the stress intensity factor 
and the experimental strengths are also in fairly good agreement, except for the 
[ ±^5 ] 2s s P ec i- mens - The deviation between the strengths predicted from the 


strains and the stress intensity factor increases with the proportion of ±45° 
pi i cs . 


CONCLUSIONS 


The method developed in reference [2] to predict fracture toughness of 
boron/aluminum laminates does not work because widespread yielding of the 
aluminum causes the compliances to be very nonlinear. An alternate method was 
developed to predict a strain intensity factor at failure. Singular strain 
fields from the Theory of Elasticity were assumed to be valid despite the 
widespread yielding. These strains are proportional to a strain intensity 
factor, just as the singular stress field is proportional to the stress 
intensity factor. A general fracture toughness parameter Q , which is 
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independent of laminate orientation, is proportional to the strain intensity 
factor at failure and another factor that depends on the elastic constants of a 
laminate. Values of strain intensity factors at failure were predicted for 
.specimens containing central crack-like slits. The specimens were made from 
five different laminate orientations: [ 0 ] 6 , [0 2 /±45) s , [±45/0 2 ) s , [ 0/±/« 5 ] g , and 

T *‘ e nonlinearity of the stress-strain curves increased with the 

proportion of ±45° plies. Failing strains were predicted from the strain 
intensity factors at failure, and strengths were predicted from the failing 
strains using uniaxial stress-strain curves. 


The predicted and measured strain intensity factors at failure and 
strengths were in generally good agreement. On the other hand, the fracture 
toughness values for laminates with 0° and off-axis plies were overestimated 
using the elastic constants. However, it was shown that, except for the [ ±45 ] 

Z o 

specimens, the fracture toughness could be predicted fairly accurately using the 
secant modulus in place of the elastic Young's modulus. This method should be 
valid for other metal matrix composites with continuous fibers. Johnson and 
Bigelow [11] using an elastic-plastic finite element model also found that 
elastic theory overestimates strength for laminates with 0° and off-axis plies. 


Although the tests and predictions were in generally good agreement, 

several discrepancies were observed: (1) The values of Q for the [0] 

c 1 6 

•specimens increased with slit length because the fiber stress at the ends of a 
slir were reduced by yielding in shear. Likewise, the experimental strengths 
were 10-20 percent greater than the predicted strengths for long slits. For 
laminates with both 0° and ±45° plies, the yielding in shear was not 
significant, and the values of Q c were independent of slit length. (2) The 

values of Q c for the [45] 2S specimens decreased with slit length. The far- 

(jeld failing strains for specimens with short slits were nearly ten times those 
for specimens with long slits. The linear-elastic analysis did not completely 
account for this difference. (3) For laminates that contained ±45° plies, the 
values of Q c increased with increasing specimen width. Strengths increased 

similarly. This effect increased with increasing proportion of ±45° plies. 
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Table 1. Tensile properties of unnotched laminates. 


Laminate 

orientation 


u 

V 

yx 

E x> 

MPa 

Hi 

b , 

yy 

NPa* 1 

n 

yy 



w - 

19.1 mm 




[°] 6 

1672.0 

0.007908 

0.2049 

143.1 

237.3 


iE3 

(V ±45 ] s 

800.1 

.007267 

.2513 


176.2 


■ 1 

[±45/0 2 ] s 

910.5 

.008205 

.2519 

134.7 

177.5 


HE 9 

[0/±45] s 

581.4 

.007008 

.2911 

129.5 

159.2 


1.743 1 

[±45] 2s 

220.6 

a . 04974 

.3247 

126.2 

126.9 


4.682 1 



W - 

50.8 mm 



( ± 45] 2s 

280.8 

a . 1511 

- 

— 



- 



w - 

101.6 mm 




MSB— 

330.9 

a . 3237 

i 

- i 

- 


- 



- Calculated with Ramberg-Osgood equation for W — 19.1 mm. 


Table II. Stress and strain intensity factors at failure. 



* 

Kq, MPaymm, for W - 

K f0 ,ymm, for 

w - 

BEBB 

50.8 mm 

101.6 mm 

19.1 mm 

50.8 mm 

101.6 mm 

t°) 6 

0.8409 

2184 

3495 

3854 

0.01091 

0.01602 

0.01794 | 

rv ±45 ) s 

. 7841 

1400 

1719 

1832 

.01251 

.01331 

.01443 

[±45/0 2 ] s 

. 7806 

' 

1699 

1821 

- 

.01348 

.01416 

f0/±45] s 

.7375 

1239 

1615 

1808 

.01454 

.01614 

.01856 

[145] 2s 

.6771 

345 

707 

703 

.02185 

.03845 

.02782 
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Figure 1.- Sketch of tensile specimen with central crack-like slit. 
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Figure 2.- Calculated stress-strain curves for the various laminate orientations. 
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Figure 3.- Calculated tangent modulus versus stress for the various laminate orientations. 
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Figure 4 - Failing strains for unnotched specimens with various laminate orientations. 
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Figure 5.- Unnotched strengths for the various laminate orientations. 
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Figure 6.- Laminate and principal lamina coordinates. 
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Figure 7.- Nondimensional parameter g versus proportion of 0°plies. 


\BAL88\HG7.BAL 






,/mm 



(b) CO2 /* 45]§ laminates. 


Figure 8.- Continued. 
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(c) [±45/0 2 ]$ laminates. 


Figure 8- Continued. 
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Figure 8 - Concluded. 
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Figure 9.- General fracture toughness parameter for the various specimen widths and laminate orientations. 
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Figure 10- Predicted and measured values of strain intensity factor at failure. 
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Figure 11- Predicted and measured values of stress intensity factor at failure. 
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Figure 12- Measured and predicted strengths for specimens with various laminate orientations. 
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(b) CO 2 /±45]g specimens. 


Figure 12.- Continued. 
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(c) it 45/O2 h specimens. 

Figure 12.- Continued. 
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Cd)[0/+4533 specimens. 


Figure 12.- Continued. 
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Figure 12.- Concluded. 
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